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Abstract 

^N I We describe quantum-octave package of functions useful for simu- 

jST ' lations of quantum algorithms and protocols. Presented package allows 

—^ , to perform simulations with mixed states. We present numerical imple- 

f^ ' mentation of important quantum mechanical operations - partial trace 

and partial transpose. Those operations are used as building blocks of 
f->i ; algorithms for analysis of entanglement and quantum error correction 

codes. Simulation of Shor's algorithm is presented as an example of 
^ ■ package capabilities. 

^- 1 Motivation 



x 



Most of software developed for simulations of quantum computing is based 
;_j ■ on finite dimensional Hilbert space formalism. |5l ITT] Mixed state model for 

^ quantum computing allows to incorporate many features which are crucial 

for analysis of entanglement and decoherence. 

Main contribution of this work is presentation of quantum-octave pack- 
age - a tool which allows to perform simulations of quantum systems using 
density operators formalism in convenient way. We present algorithms used 
in implementation of partial operations in quantum-octave. Sample results 
obtained using quantum-octave present capabilities of package. 
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2 Implementation in GNU Octave language 

Package quantum-octave is implemented in GNU Octave language. [3] List 
of functions and their detailed description can be found on project web 
page.fTU] We present here only the most important features of quantum-octave. 

Most of quantum-octave functions can be applied to density operators, 
which are represented by normalise matrices. Package offers also functions 
for operations on pure states. Base ket vectors and their linear combi- 
nations are easily constructed and transformed into corresponding den- 
sity operators using State function, and it is easy to obtain mixtures of 
states using MixStates function. Computation process is performed using 
Evolve (gate, state) function, where gate is unitary operation. 

Lowlevel functions BinVec2Dec and Dec2BinVec allow for convenient use 
of register notation (|1) (gi |0) |0) . . . (gi |0)). They are implemented in C++ 
using liboctave are they are used for manipulations of registers in PTrace 
and PTranspose functions. For example sequence of operaions 

y = Dec2BinVec (3,4) ; 

tmp = y(2) ; 

y(2) = y(3); 

y(3) = tmp; 

X = Ket(BinVec2Dec(y)) ; 

allows to permute second and third qubits in base state in four dimensional 
space. 

3 Algorithms for partial operations 

Partial operations are crucial for analysis of quantum information process- 
ing. Partial transposition allows to distinguish entangled and sparable states 
in low dimensional systems. 4 This operation was used in Ref. J2| for con- 
structing computable measure of entanglement - negativity.^ 

3.1 Partial transposition using registers 

In quantum-octave partial transposition is based on permutation of qubits 
indexes. Suppose we perform partial transposition with respect to some 



^Negativity is defined as a sum of negative eigenvalues of density matrix after partial 
transposition. It is easy to clieck tfiat tliis definition is equivalent to formula Af{p) — 

■Li^ — j ' '- ~ , where ||p||i denotes trace norm of density matrix p. 



qubits in input state. Let Q = {gi, q2, ■ ■ ■ , Qn} C N be the set of numbers 
labelling those qubits. Using register notation one can write 

{a\p^n'-'^- 1/3) = {Via, /3, gi, . . . , qn)\p\Vif3, a,qi,..., Qn)), (1) 

where 

and ai, l3i,Vi represent i-th number in binary representation of a, /5 and V. 

3.2 Partial trace 

For density matrix p £ SiTiA^'HB) one can obtain description of subsytems 
A and B using partial trace operation. [^1 15] In quantum- octave this oper- 
ation is implemented as follows. Let Q = {qi, q2, ■ ■ ■ , qn} be a given set of 
number labelling qubits. Density operator describing state of this qubits is 
defined as matrix with elements 

1 1 

(a|trgi,,,,...,,„(p)|/3) = ^ ... Y, i^i^'ku . . . ,kn)\p\WiP,h, . . . ,kn)), 

ki=0 kn=0 

(3) 
where a and /3 are binary digits of length m = N — n. Number ki, . . . ,kn G 
{0, 1} label base vectors in Hilbert spaces of qubits qi, . . . ,qn respectively. 
Numbers W are constructed as follows 

W{a,ku...,kn) = {WuW2,...,WN),W, = (''J^ '.^^ . (4) 

In quantum-octave construction of numbers W is performed using func- 
tion BuildBinaryVector implemented in C++. This is motivated by ex- 
tensive usage of this function. If one wants to perform partial trace oper- 
ation on A^-qubit system with respect to its m-qubit subsystem, function 
BuildBinaryVector have to be called 2^(^~™)^ times. 

One should note that this algorithm allows to perform partial trace with 
respect to any set of qubits. 

4 Shor's algorithm with mixed sates — fideUty and 
distances measures 

In this section we present sample results obtained by using quantum-octave. 
We analyse influence of states mixing on Shor's algorithm.^ The idea was 
^The subject of analysis is restricted to quantum part of Shor's algorithm. 
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(a) Fidelity between density matrices as a func- 
tions of 'purity' of state 
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(b) Fidelity between probability distributions as 
a functions of 'purity' of state 
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Figure 1: Comparison odistancece measures between density matrices 
(Fig. 1(a)) and between probability distributions (Figs. l(b)|l(d) ) for sim- 
ulation of Shor's algorithm 



to observe how quantum state, obtained by performing quantum algorithm 
on non-pure state, differs from state obtained in ideal (pure state) case. On 
each step of simulation Initial state pp was prepared as mixture of pure state 
po = I 00 . . . 0) (00 ... I and maximally mixed state Id/d 
d d 



Pp = {I -p)po+pId/d, 



(5) 



where p G {0.00, 0.01, 0.02, . . . , 1.00} and Id is d-dimensional identity matrix. 
For each state obtained in such way the trivial and Shor's circuits were ap- 
phed and for each obtained outcome two distances and two fideUty measures 
were calculated. Results of those simulations are presented in Fig. 1. Shor's 
algorithm is performed on 7 qubits, but only the first 3 qubits are taken 
into account during calculations of fidelity and distance measures because 
the last 4 qubits decohere and normally aren't measured. 

Fidelity between density matrices, as given in Ref. IHl is proper mea- 
sure which shows how much two mixed states differ. jHI To compare different 
methods of calculating distance between states we have chosen three func- 
tions (Eqs. 13 IHl and EJ that operate on probability distributions. Proba- 
bility distributions are obtained by performing measurement of observable 
Z ® Z ® Z. Formulas used to calculate measures are presented below. Let 
Pi,P2 be probability distributions, and pi,P2 ~ density matrices. 

• Fidelity beetween density matrices 



F{pi,p2)=tr^^P2^/p^, (6) 

Fidelity beetwen probability distributions 

F{pi,P2) = ^ ^/Pi{x)^/p2{x), (7) 



xex 



where p = \{pi + P2)'-, 

2( . V^ {pi{x)-p{x)f 

• Trace-distance 

F{pi,P2) = ^^\pi{x)-p2(.x)\. (9) 

One can conclude that for small addition of noise results differs wery much 
from the ideal case. The loss of quality grows slower for bigger contribution 
of noise. 



5 Final remarks 

Package quantum-octave allows to perform simulations of mixed states 
quantum computation in convenient way. It also provides functions for anal- 
ysis of entanglement and quantum errors. Quantum protocols and algorithm 
such as teleportation,^ qubit authentication!^ or Shor's algorithml^I can 
be easily implemented as GNU Octave functions and the use in futher simu- 
lations. This allows to perform simulations operating on high level quantum 
primitives (e.g. qubit authentication, error correcting code). 
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